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Abstract 

We present results of a direct search for the decay Ks — > e^e~ with the KLOE de- 
tector, obtained with a sample of e~^e~ ^ </> — > KsKl events produced at DA<I>NE, 
the Frascati (/)-factory, for an integrated luminosity of 1.3 fb~^. The Standard Model 
prediction for this decay is BIi{Ks e~^e~) = 1.6 x 10~^^. The search has been 
performed tagging the Ks decays by simultaneous detection of a interaction 
in the calorimeter. Background rejection has been optimized by using both kine- 
matic cuts and particle identification. At the end of the analysis chain we find 
BR(ir5 e+e~) < 2.1 x 10"® at 90% CL, which improves by a factor of ~ 7 on 
the previous best result, obtained by CPLEAR experiment. 



1 Introduction 



The decay Ks — >■ e'*"e~, like the decay Kl — ^ e+e^ or Kl yW^/W , is a 
flavour-changing neutral-current process, suppressed in the Standard Model 
and dominated by the two-photon intermediate state [1]. For both Ks andKi, 
the e~^e~ channel is much more suppressed than the fi^fi^ one (by a factor of 
~ 250) because of the e — /j, mass difference. The diagram corresponding to 
the process Ks — > 7*7* — ^ is shown in Fig. 1. Using Chiral Perturbation 




Fig. 1. Long distance contribution to Ks i^i process, mediated by two-photon 
exchange. 

Theory (xpT) to order O(p^), G. Ecker and A. Pich evaluated the ratio [1] 
T{Ks e+e-)/T{Ks ^ 77) = 8 x 10"^ with 10% uncertainty. Using the 
present average [2] for BIi{Ks — ^ 77) we obtain the Standard Model prediction 
BIi{Ks — > e+e~) ~ 10~^^. A value significantly higher than expected would 
point to new physics. The best experimental limit for BR{Ks e^e~) has 
been measured by CPLEAR [3], and it is equal to 1.4 x 10"^ at 90% CL. 



2 



Here we present a new measurement of this channel, which improves on the 
previous result by a factor of ~ 7. This paper is organized as follows: in the 
next section, a brief description of the KLOE experimental setup is given; 
in section 3 the selection criteria for the decays of interest are summarized. 
Results are presented in section 4. 



2 Experimental setup 



The data were collected with KLOE detector at DA$NE, the Frascati (p- 
factory. DA$NE is an e^e~ collider that operates at a center-of-mass energy 
of ~ 1020 MeV, the mass of the meson. Positron and electron beams of 
equal energy collide at an angle of vr — 25 mrad, producing mesons with a 
small momentum in the horizontal plane: ~ 13MeV/c. mesons decay 
~ 34% of the time into nearly coUinear K^K^ pairs. Because J^'"{(f)) = 1 , 
the kaon pair is in an antisymmetric state, so that the final state is always 
KsKl. The contamination from K^Ki and KgKs final states is neglegible. 
Therefore, the detection of a Kl signals the presence of a Ks of known mo- 
mentum and direction, independently of its decay mode. This technique is 
called Ks tagging. The analyzed sample corresponds to an integrated lumi- 
nosity of ~ 1.3 fb^^, yielding ~ 1.4 billion of KsK^ pairs. The KLOE detector 

m ^ m n m m * 

II ' '-^ ' . YOKE > ^ / ' ' I ^ 
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Fig. 2. Vertical cross section of the KLOE detector. 

(Fig. 2) consists of a large cylindrical drift chamber (DC), surrounded by 
a lead/scintillating- fiber sampling calorimeter (EMC). A superconductig coil 
surrounding the calorimeter provides a 0.52 T magnetic field. The drift cham- 
ber [4], which is 4 m in diameter and 3.3 m long, has 12582 all-stereo tungsten 
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sense wires and 37746 aluminium field wires. The chamber shell is made of 
carbon-fiber/epoxy composite, and the gas used is a 90% helium, 10% isobu- 
tane mixture. These features maximize trasparency to photons and reduce 
Kl^Ks regeneration and multiple scattering. The DC position resolutions are 
(Txy ~ 150//m and (Tz ~ 2 mm. DC momentum resolution is a{p±)/p± ~ 4%. 
Vertices are reconstructed with a spatial resolution of ~ 3 mm. 

The calorimeter [5] is divided into a barrel and two endcaps, contains a total 
of 88 modules, and covers 98% of the solid angle. The modules are read out 
at both ends by photomultiplier tubes. The arrival times of particles and the 
three-dimensional positions of the energy deposits are determined from the 
signals at the two ends. The read-out granularity is ~ 4.4 x 4.4 cm^; fired 
"cells" close in space and time are arranged into a "calorimeter cluster". For 
each cluster, the energy Ed is the sum of the cell energies, and the time 
tci and the position are calculated as energy-weighted averages over the 



fired cells. The energy and time resolutions are (Je/E = 5.7%/ J E(GeV) and 



The calorimeter trigger [6] requires two local energy deposits above a threshold 
of 50 MeV in the barrel and 150 MeV in the endcaps. Recognition and rejec- 
tion of cosmic-ray events is also performed at the trigger level: events with 
two energy deposits above a 30 MeV threshold in the outermost calorimeter 
plane are rejected as cosmic-ray events. Moreover, to reject residual cosmic 
rays and machine background events, an offline software filter (FILFO) ex- 
ploits calorimeter and DC informations before tracks are reconstructed [7]. 
The trigger has a large time spread with respect to the beam crossing time. 
However, it is synchronized with the machine RF divided by 4, Tgync ~ 10.8 ns, 
with an accuracy of 50 ps. An estimate of the time of the bunch crossing pro- 
ducing an event is determined offline during event recostruction. This value 
is subtracted from the measured cluster times to obtain particle time-of-fflght 
(TOF) measurements. 

The response of the detector to the decays of interest and the various back- 
grounds were studied by using the KLOE Monte Carlo (MC) simulation 
program [7]. Changes in the machine operation and background conditions 
are simulated on a run-by-run basis. The most important parameters are the 
beam energies and the crossing angle, which are obtained from the analysis of 
Bhabha scattering events with polar angle above 45 degrees. The average 
value of the center-of-mass energy is evaluated with a precision of 30keV for 
each 100 nb~^ of integrated luminosity. To study the background rejection, a 
MC sample of decays to all possible final states has been used, equivalent to 
an integrated luminosity of ~ 1.3 fb~^. A MC sample of ~ 45000 signal events 
has been also produced, to measure the analysis efficiency. 




at = 57 ps/ JE{GeV) © 100 ps, respectively. 
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3 Data analysis 



3.1 Ks tagging 

The identification of fr^-interaction in the EMC is used to tag the presence of 
Ks mesons. The mean decay lenghts of Kg and are A5 ~ 0.6 cm and \l ~ 
350 cm, respectively. About 50% of Kls therefore reach the calorimeter before 
decaying. The Kl interaction in the calorimeter barrel (i^^crash) is identified 
by requiring a cluster of energy greater than 125 MeV not associated with any 
track, and whose time corresponds to a velocity j3 = rd/ctd compatible with 
the kaon velocity in the center of mass, /?* ~ 0.216, after the residual (p 
motion is considered. Cutting at 0.17 < (3* < 0.28 we selected ~ 450 miUion 
Xg-tagged events (i^^crash events in the following), which are used as a starting 
sample for the Ks e^e~ search. 



3.2 Signal preselection and background normalization 

Ks e~^e~ events are selected by requiring the presence of two tracks of 
opposite charge with their point of closest approach to the origin inside a 
cylinder 4 cm in radius and 10 cm in length along the beam line. Moreover, 
the two tracks are required to form a vertex with position in the transverse 
plane p < 4 cm. The track momenta and polar angles must satisfy the fiducial 
cuts 120MeV/c <p< 350MeV/c and 30° < ^ < 150°. The tracks must also 
reach the EMC without spiralling, and have an associated cluster. In Fig. 3, the 
two-track invariant mass evaluated in electron hypothesis (Mge) is shown for 
both MC signal and background samples. A preselection cut requiring Mg^ > 
420MeV/c^ has been apphed, which rejects most of Ks — > n'^n' events, for 
which Mee ~ 409 MeV/c^. The residual background has two main components: 
Ks Ti^Ti~ events, populating the low M^e region, and — > 7r"'"7r~7r° events, 
spreading over the whole spectrum. The Ks tt+tt" events have such a wrong 
reconstructed M^e because of track momentum resolution or one pion decaying 
into a muon. The — > 7r+7r~7r° events enter the preselection because of a 
machine background cluster, accidentally satisfying the -f^crash algorithm. After 
preselection we are left with ~ 10^ events. To have a better separation between 
signal and background, a x^-like variable is defined, collecting informations 
from the clusters associated to the candidate electron tracks. Using the MC 
signal events we built likelihood functions based on: 

• the sum and the difference of 5t for the two tracks, where 5t = td — L / I3c is 
evaluated in electron hypothesis; 

• the ratio E/p between the cluster energy and the track momentum, for both 
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Fig. 3. Two-track invariant mass evaluated in electron hypothesis for MC signal 
(red) and background (black) events. 

charges; 

• the cluster centroid, for both charges. 

In Fig. 4, the scatter plot of versus Mge is shown, for MC signal and 
background events. The spectrum for background is concentrated at higher 
values respect to signal, since both Ks tt+tt^ and (p tt^tt^tt^ events have 
pions in the final state. To assess the MC background normalization, two 
sidebands are defined in the invariant mass: Mee < 460MeV/c^ (region 1 in 
Fig. 4), and M^e > 530 MeV/c^ (region 3 in Fig. 4). Ks vr+vr" events largerly 
dominate on ^ tt+tt^tt^ in region 1, the opposite occurring in region 3. In 
the two sidebands, a normalization factor is evaluated for each background 
component, by fitting the MC spectra to data. A comparison between data 
and MC after the fit is shown in Fig. 5. 



3.3 Background rejection 

A signal box to select the Ks e^e~ events can be conveniently defined in 
the region 2 of the M^e — plane (see Fig. 4); nevertheless we investigated 
some more independent requirements in order to reduce the background con- 
tamination as much as possible before applying the Mee ~ selection. These 
cuts have been tuned on the Mee sidebands, which are also used to check 
data-MC consistency after each step of the analysis, and they are summarized 
below. Charged pions from Ks —>■ tt^tt^ decay have a momentum in the Ks 
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Fig. 4. Scatter plot of versus M^e for MC signal (red) and background events 
(black) 
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Fig. 5. Data-MC comparison for M^e spectra in region 1 (left) and region 3 (right), 
after normalization: data are represented by the black solid line, MC by red points. 

rest frame p* ~ 206MeV/c. The distribution of track momenta in the Ks 
rest frame, evaluated in the pion mass hypothesis, is shown in Fig. 6, for MC 
background in region 1, and for MC signaL For most of Ks —>■ vr+vr" decays, 
at least one pion has well reconstructed momentum, so that the requirement 
min{p*{l),p*{2)) > 220MeV/c rejects ~ 99.8% of these events, while retain- 
ing ~ 97% of the signal. To reduce the 7r+7r~7r° background, we can 
use the fact these events have two photons coming from the interaction point. 
The distribution of the number of prompt photons {St = — Vd/c < Sa^) 
for data and MC background in region 3, and for MC signal, is shown in 
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Fig. 7. We require Np^ompt < 1, thus rejecting ~ 65% of background and only 
~ 0.1% of signal events. Further rejection on ^ tt^tt^tt'^ events is achieved 
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Fig. 7. Number of prompt photons (Nprompt) for data (black) and MC background 
(red) in region 3 (left), and for MC signal (right). 



by cutting on the missing mass, evaluated as M^^^^ = (^P^ — P+ — V-j where 

is the four-momentum and p± are the charged track four-momenta, in 
pion hypothesis. Distribution of M^iss is shown in Fig. 8, for data and MC 
background in region 3, and for MC signal. We cut at M^iss > 380MeV/c^, 
rejecting almost completely the 3-pion background, peaking at the 7r° mass. 
A comparison between data and MC counts on the Mee sidebands after each 
cut is shown in Tab. 1, which demonstrates the reliability of the background 
simulation. 
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Fig. 8. Missing mass {Mmiss) for data (black) and MC background (red) in region 
3 (left), and for MC signal (right). 



Cut 


Region 1 
Data MC bkg A/a 


Region 3 
Data MC bkg A/cr 


* 

-^prompt 


6061 6285(125) -1.5 
2778 2982(80) -2.1 
1370 1407(53) -0.6 


13144 13590(184) -2.0 
4672 4745(103) -0.6 
20 5(3) 2.8 



Table 1 

Data and MC counts in the Mee sidebands after each step of the background rejec- 
tion; the relative difference between data and MC is also reported. 



3.4 Signal box definition 



A signal box is defined in M^e — plane. To achieve the best background 
rejection, with maximum signal efficiency, we scan over a large set of Mee — 
cut configurations. This optimization procedure is based on MC only, without 
looking at the number of observed events (Nobs) on data. The best configura- 
tion is: 

492 MeV < Mee < 504 MeV 



Applying this selection to the data sample we obtain Nobs = 3. The expected 
background estimated from MC is fiB = 7.1±3.6, which takes into account the 
MC statistics and uncertainty on the normalization factors. Using a bayesian 
approach [8], we evaluate the upper limit on the expected number of signal 
events (Xs to be UL{fis) = 4.3, at 90% CL. 
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3.5 Radiative corrections 



Given the chosen invariant mass selection, we actually measure the upper 
limit on Ks — > e+e~(7), with E* < 6MeV. Two processes are expected to 
contribute to photon emission, not interf erring with each other: 

• the inner bremsstrahlung photon emission, Ks e^e~ + 77^; 

• a Kg — > 77 decay, with one photon conversion, Ks — > 77* — > ^e^e~ . 

For the first process we estimate a correction given by [9]: 

r(X5-e+e-(7)) 

The second process is strongly peaked around Mee ~ Im^ [10], and a BR{Ks 
77* 7e+e~) ~ 2 X 10~^^ is expected in our final mass window. So, our result 
consists in a limit on the IB emission, being insensitive, in this range, to 
the photon conversion process. 



4 Results 



The total selection efficiency on Ks — > e^e events is evaluated by MC, using 
the following parametrization: 

where e{Kcrash) is the tagging efficiency, and €{sele\K crash) is the signal selec- 
tion efficiency on the sample of tagged events. The number of Ks 7r~^TT~ 
events N^+t^- counted on the same sample of Ks tagged events is used as 
normalization, with a similar expression for the efficiency. The upper limit on 
BR{Ks — > e+e~) is evaluated as follows: 



UL{BR{Ks e+e-)) = 



jjj, X ^ e^+^-{sele\K crash) BR{Ks ^ n'^ir' 

U L{Hs) X I'Ctag X — — — X 



^sig^^^^^\K crash) -^i 



where TZtag is the tagging efficiency ratio, corresponding to a small correction 
due to the i^^crash algorithm dependence on Ks decay mode. Determination of 
both TZtag and e^+^-{sele\Kcrash) is discussed in detail in Ref. [11]. Using as 
input values: 

• UL{i^s) = 4.3, at 90% CL. 
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• esig{sele\Kcrash) ^ Om? x 0^ = 0.558(4) 

cuts rad 

• Utag = 0.9634(1) 

• e^+^-{sele\K crash) = 0.6102(5) 

• N^+^- = 148, 184, 688 

we obtain 

UL{BR{Ks e+e-(7))) = 2.1 x 10-^ at 90% CL. 

Systematics uncertanties, related to background normalization, are at the level 
of 2%. Our measurement improves by a factor of ~ 7 on the CPLEAR re- 
sult [3], for the first time including radiative corrections in the evaluation of 
the upper limit. 



References 



[I] G. Ecker and A. Pich, Nucl. Phys. B 366 (1991) 189. 
[2] W.-M. Yao et al., J. Phys. J. G 33 (2006) 1. 

[3] A. Angelopoulos et al., Phys. Lett. B 413 (1997) 232. 

[4] KLOE collaboration, M. Adinolfi et al, Nucl. Istrum. Meth. A 488 (2002) 51. 

[5] KLOE collaboration, M. Admolfi et al, Nucl. Istrum. Meth. A 482 (2002) 363. 

[6] KLOE collaboration, M. Adinolfi et al, Nucl. Istrum. Meth. A 492 (2002) 134. 

[7] KLOE collaboration, F. Ambrosino et al, Nucl. Istrum. Meth. A 534 (2004) 
403. 

[8] O. Helene, Nucl. lustrum. Meth. 212 (1983) 319. 
[9] C. Gatti, Eur. Phys. J. C 45 (2005) 417. 

[10] L. Bergstrom, E. Masso and P. Singer, Phys. Lett. B 131 (1983) 229. 

[II] KLOE coUaboration, F. Ambrosino et al, Eur. Phys. J. G 48 (2006) 767 



11 



